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Abstract 

We study the baryonic charmonium decays of B mesons, B + — > rj c K + and B + — > J/ipK + , 
where the r\ c and J/ip subsequently decay into a pp or AA pair. We measure the J/tp — > pp and AA 
anisotropy parameters, as = —0.60 ± 0.13 ± 0.14 (pp), —0.44 ± 0.51 ± 0.31 (AA) and compare to 
results from e + e~ — > J ftp formation experiments. We also report the first observation of r/ c — ► AA. 
The measured branching fraction is B(i] c -» AA) = (0.87lg;|f(5*at)l8;S|(sj/s*)±0.27(PDG)) x 10~ 3 . 
This study is based on a 357 fb _1 data sample recorded on the T(4S) resonance with the Belle 
detector at the KEKB asymmetric-energy e + e~ collider. 

PACS numbers: 13.25.Gv, 14.40. Gx, 13.40.Hq 
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There have been many observations of baryonic three-body B decays in recent years 0, 0, 
IH 0, H[ . An interesting feature of these observations is the presence of peaks near threshold 
in the mass spectra of the baryon-antibaryon pair. These enhancements are not likely to be 
resonance states, as the baryon angular distributions are not symmetric in their respective 
helicity frames [6j]. Other visible structures in the mass spectra arise from charmonium 
decays. It is natural to compare the baryon angular distributions from charmonium decays 
with those in the region of the threshold enhancement. There is particular interest in 
J/ijj — > pv, where the proton an gula r distribution has been studied by DASP [7] , DM2 , 
Mark I Q, Mark II @ and BES lUHIljj. W 

mesons from the reaction e + e — > J/ip are 
produced predominantly in helicity ±1 states. Therefore, the baryon angular distributions 
are proportional to 1 + a cos 2 6, where 9 is the baryon polar angle in the J/ip helicity frame. 
Many theoretical predictions |l4| exist for the value of a. 

Study of two-body baryonic decays of charmonia at a 5-factory has several different 
features as compared with an e + e~ machine running at the J ftp mass. J/ip mesons from the 
two body decay of B mesons accompanied by spin zero particles are in a pure helicity zero 
state. This provides a useful cross check for previous measurements. The charmonia from 
B decays do not suffer from poor acceptance near the beam pipe, and events with | cos#| 
near 1 can be detected. Such events are very effective for determining a. Requiring that 
the J/ip originate from a B decay eliminates e + e~ — > qq — > pp background, where q stands 
for auorrf quark. For e + e~ — > J/ip — > PP, this background cannot be separated from the 
signal on an event-by-event basis. 

In the study of two-body baryonic decays of charmonia we focus on the decay processes 
B + — > ppK + and B + — > AAK + ^5[. We report the first observation of r\ c — > AA. There is 
little information about r] c decays into baryon-antibaryon pairs except for i] c — > pp. Measur- 
ing decay rates of the rj c to different baryon-antibaryon modes is a useful check for theoretical 
predictions 0] and can shed light on quark-diquark dynamics. 

We use a 357 fb _1 data sample consisting of 386 x 10 6 BB pairs collected by the Belle 
detector at the KEKB asymmetric energy e+e" (3.5 on 8 GeV) collider [17]. The Belle 
detector is a large solid angle magnetic spectrometer that consists of a four layer silicon 
vertex detector (SVD), a 50 layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time of flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of Csl (Tl) crystals located inside 
a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux return 
located outside of the coil is instrumented to detect K\ mesons and to identify muons. The 
detector is described in detail elsewhere 

The event selection criteria are based on information obtained from the tracking system 
(SVD+CDC) and the hadron identification system (CDC+ACC+TOF). We follow the same 
procedure as in ref. 0] to select proton and kaon candidates. A candidates are reconstructed 
via the pir~ channel using the method described in ref. [lit ]. 

To identify the reconstructed B meson candidates, we use the beam energy constrained 



mass Mb c = y E^ eam — p\ and the energy difference AE = E B — E heiim , where E hes , m is the 
beam energy, and pb and Eb are the momentum and energy of the reconstructed B meson in 
the rest frame of the T(4S). The signal region is defined as 5.2 GeV/c 2 < Mb c < 5.29 GeV/c 2 
and -0.1 GeV < AE < 0.2 GeV. The signal peaks at M bc = 5.279 GeV/c 2 and AE = 0. 

The dominant background arises from continuum e + e~ — > qq processes. The background 
from b —>■ c and from B decays into charmless final states is negligible. In the T(4S) rest 
frame, continuum events are jet-like while BB events are more spherical. The reconstructed 
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momenta of final state particles are used to form several event shape variables (e.g. thrust 
angle, Fox- Wolfram moments, etc.) in order to categorize each event. We follow the scheme 
described in ref. |2o| that combines seven event shape variables into a Fisher discriminant 
to suppress continuum background. 

Probability density functions (PDFs) for the Fisher discriminant and the cosine of the 
angle between the B flight direction and the beam direction in the T(4S) rest frame are 
combined to form the signal (background) likelihood £ s (&). The signal PDFs are determined 
from GEANT based Monte Carlo (MC) simulation and the background PDFs are obtained 
from sideband data with My, c < 5.26 GeV/c 2 . We require the likelihood ratio 1Z = L s j (C s + 
Cb) to be greater than 0.4 for both ppK + and AAf^ + modes. These selection criteria suppress 
approximately 69% (66%) of the background while retaining 92% (91%) of the signal for the 
ppK + (AAK + ) mode. If there are multiple B candidates in an event, we select the one with 
the best y 2 value from the B decay vertex fit. Multiple B candidates are found in less than 
2% (5%) of events for the ppK + (AAK + ) mode. 




M - (GeV/c 2 ) M M (GeV/c 2 ) 



FIG. 1: (a) B signal yield versus M p p and (b) B signal yield versus M AA . The inset shows the 
r] c -J/ip mass region. The curves represent the unbinned likelihood fits to the data. 



We use an unbinned extended maximum likelihood fit to estimate the B signal yield. 
For the signal PDF, we use a Gaussian in M^ c and a double Gaussian in AE. We fix the 
parameters of these functions to the values determined from MC simulation jilj . Background 
shapes are fixed from fitting to the sideband events in the region 3.14 GeV/c 2 < M p p < 3.34 
GeV/c 2 . The Mt, c background is modeled using a parametrization used by the ARGUS 
collaboration 2^| . The AE background shape is modeled by a first order polynomial. 

We determine B signal yields in 10 MeV/c 2 wide M p p (Maa) m ass bins from the kinematic 
threshold to 4.5 GeV/c 2 ; the result is shown in Fig. da) (Fig. ^b)). There are clear r\ c and 
J/i/j peaks in the mass spectrum. We use a relativistic Breit-Wigner function for the r] c peak, 
a Gaussian for the J/ip peak, and a linear function for the non-resonant background. The 
Breit-Wigner function is convolved with the detector response function, which is taken from 
the J/ip peak. A maximum likelihood fit to the data is shown in the inset. We obtain an 
r] c mass of M^ = 2971 ± 3±j MeV/c 2 (2974 ± 7t? MeV/c 2 ) and a width of F(r) c ) = 48±? ± 5 
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MeV/c 2 (40 ± 19 ± 5 MeV/c 2 ) for the r] c — > pp (i] c — > AA) mode. The systematic errors are 
determined from the differences of J /if) peaks between data and the PDG |25[ value, and 
by varying different fit shapes for r] c signal and background assuming no interference effect 
between them. The width is larger than the PDG average but is consistent with recent 
BaBar 0, |23| and previous Belle j24[ measurements. 
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FIG. 2: (a) Likelihood fit and (b) x 2 fit results of the J /if) — > pp helicity angle distribution. In 
the maximum likelihood fit plot, the solid, dotted solid, and dashed line represent the combined 
fit result, fitted signal, and fitted background, respectively. In the x 2 fit plot, the inset shows the 
fit result for B signal yield of B + — > J/tfjK + , J /if) — > ■ 

We define the J /if) signal region as 3.075 GeV/c 2 < M pP < 3.117 GeV/c 2 and use events 
in this signal region to study the proton angular distribution in the helicity frame of the 
J/if>. The helicity angle 9x is defined as the angle between the proton flight direction and the 
direction opposite to the flight of the kaon in the J / ip rest frame. The angular distribution 
of J /if) in the helicity zero state is parameterized as P(aB,cos9x) = (1 + «b cos 2 9x)/(2 + 
2as/3) with as = — 2a/(a+l) 26J. Here, a is the anisotropy parameter determined 
from the angular distribution of J/if> in helicity= ±1 states produced in e + e~ — > J/if>- 
Previous measurements 0, 0, H, Efl EH 03 give an average of a = 0.66 ± 0.05 for pp 
and a = 0.65 ± 0.11 for AA; these values correspond to as = —0.80 ± 0.04 for pp and 
a B = -0.79 ± 0.08 for AA. 

For analysis of the angular distribution, we define a likelihood L, 

e -(N 3 +N b ) N 

— H[N s P s (M hCi , AEi)e(cos6xi)P(a B , cosB Xi ) 

i=l 

+N b P b (M hcv AEi,cos6 Xi )}, 

where ag is a fit parameter in addition to and N b , e(cos9 x ) is the efficiency function and 
eP is normalized to 1. The efficiency e(cos8x) obtained from the signal MC simulation is 
flat as a function of cos#x- From a study of signal MC simulation, we find that there is no 
correlation between Mb c , AE and 6x- The background PDF as a function of Mt, c , AE, and 
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cos 6x is determined from M pP sideband data. Fig. 12(a) shows the result of the fit to the 
J/ijj — > pp candidates in the entire Mt, c , AE region. We determine as to be — 0.60±0.13. As 
a cross-check, we use a \ 2 method and fit the efficiency corrected B signal yields in bins of 
cos 9x to a 1 + a>B cos 2 9x parametrization. The results of the fit are shown in Fig. E(b). We 
obtain «b = —0.53 ± 0.15 with \ 2 /d.o.f. = 0.9, consistent with the result of the unbinned 
fit. We measure the angular distribution of J/ip — > ji + fi~ decays from B + — > J/ipK + to 
verify the fitting procedure. The result is shown in the inset of Fig. 12(b) . The fitted value 
agrees with the expectation for massless fermions. 

We determine the systematic error in by varying the value of various selection cuts and 
parameters of PDFs to check for trends in the value of as- These trends are parametrized 
by a linear function. We then quote the change in «b along the line between the selected 
point and the far end of the tested region as a systematic error. Note that this is a conser- 
vative estimate, since statistical fluctuations also contribute to changes in as- We assign 
a systematic error of 0.08 for the 1Z selection, 0.06 for proton/kaon selection, and 0.02 for 
fitting PDFs. Other systematic errors are negligible. The observed difference between the 
maximum likelihood method and the x 2 method is also included in the systematic error. 
The total systematic uncertainty in as is 0.13. 

There are several complicating factors in the analysis of B + — > AAK + decays relative to 
B + — > ppK + decays. The efficiency for detecting slow pion from A decays is small. As a 
result, the A reconstruction efficiency is non-uniform as a function of polar angle (cos^p) of 
the secondary proton in the A helicity frame, and is correlated with cos9x, where X refers 
to the A. The likelihood function is similar to the previous one except that the angular part 
contains two more variables, cos 9 P and cos 9 p . The efficiency function e(cos 9x, cos 9 P , cos 9p) 
is obtained from a signal MC sample with 4 x 10 6 events. The background PDF is determined 
from M AA sideband data in the region 3.14 GeV/c 2 < M AA < 3.54 GeV/c 2 . The value of 
as obtained from the maximum likelihood fit is —0.44 ± 0.51 ± 0.31, where the systematic 
error is determined from the same procedure as that used for J/ip — > pp decays. 

We define an rj c signal region as 2.94 GeV/c 2 < M AA < 3.02 GeV/c 2 . Signal peaks are 
visible in the M^ c and AE distributions. The fitted B signal yield, efficiency, and obtained 
branching fraction are shown in Table HJ The maximum likelihood fit for B + — > i] c K + ,i] c — > 
AA gives a yield of 19.5l|;4 with a statistical significance of 7.9 standard deviations. The 

significance is defined as y— 21n(Lo / L mffl ) , where Lq and L max are the likelihood values 
returned by the fit with the signal yield fixed to zero and at its best fit value, respectively. 
The fit yield is consistent with the yield (18.2 ± 4.8) obtained from the first fit shown in 
Fig. IHb). As a cross check, the obtained £>(J/V —>■ pp,AA) are in good agreement with the 
world average and with the latest BES result [l3j . We also determine the branching fraction 
ratios: B(r) c -> AA)/B(r] c -> pp) = 0.67±8;i§ ± 0.12 and B{J/i/; -> AA)/B(J/V> -> pp) = 
0.90i ; 14 ±0.10, where common systematic errors in the numerator and denominator cancel. 
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TABLE I: Measured branching fractions for J/ip,r] c —> pp, AA 



Modes 


Yield 


Eff.(%) 


B.F.Product(l(T 6 ) 


B(J/iP, Vc ^pp,AA)(10- 3 ) 


B + — > f] c K + ,r] c — > pp 
B+ -> ?? C K + , r? c -> AA 
B+ -> J/^+, J/^ -> 

fi+ -> J/v> -> aA 


195+11 
19.5133 

3171H 
45.9±£? 


35.8lgJ 

5.312:1 

37.3t°j 
c q+o.3 

°- y -0.3 


-i ^9+0.11+0.16 

i-^-O.ll-O^O 
n QC-+0.25+0.08 

u.yo_ 2 2_o n 

2.2i±8:il =b 6.10 
2.00±g;|| ± 0.34 


1.58 ± 0.12+^ ± 0.47 a 
0.871^1^ ±0.27 b 
2.21 ± 0.13 ± 0.31 ±0.10 c 
2.00±g;H ± 0.34 ± 0.08 c 



a -> r/ c X+) = 0.9 ± 0.27 x lO" 3 ^. 

b We use -> r) c K+,ri c -> kk)/B(B + rj c K+,r) c -> pp) = 0.671^6 ± 0.12 measured in this paper 

and B(?7 C -> pp) = 1.3 ± 0.4 x 10" 3 0. 
c B(B+ -> J/i)K+) = 1.00 ± 0.04 x lO- 3 0. 



Systematic uncertainties are studied using high statistics control samples. For proton 
identification, we use a A — > p-n~ sample, while for K/ir identification we use a D* + — > D°7r + , 
Z}° — > K~n + sample. The tracking efficiency is studied with fully and partially reconstructed 
.D* samples. The modeling of the TZ continuum suppression requirement is studied with the a 
topologically similar control sample, B + — > J/ipK + , J/ip — > For A reconstruction, we 

have an additional uncertainty on the efficiency for detecting tracks away from the IP. The 
size of this uncertainty is determined from the difference between A decay-time distributions 
in data and MC simulation. Based on these studies, we assign a 1% error for each track, 2% 
for each proton identification, 1% for each kaon/pion identification, an additional 3% for A 
reconstruction and 3% for the 1Z selection. 

The systematic uncertainty in the fit yield is studied by varying the parameters of the 
signal and background PDFs and is approximately 5%. The MC statistical uncertainty and 
modeling contributes a 5% error. The error on the number of BB pairs is determined to be 
1%, where the branching fractions of T(4S) to neutral and charged BB pairs are assumed to 
be equal. The non-charmonium feed-down background below the r] c mass region is estimated 
to be 8% and 12% for the pp and AA modes, respectively. 

The correlated errors are added linearly and combined quadratically with the uncorrelated 
errors in the systematic error calculation. The total systematic uncertainties are 14% and 
17% for the ppK + , and AA.fr + modes, respectively. 

In summary, using 386 x 10 6 BB events, we measure the branching fractions of J/ip — > pp, 
r] c — > pp, J ftp — > AA and rj c — > AA from B + — > ppK + and B + — > AAK + decays. We measure 
the parameter as for baryonic J/ip decays. The parameters oib are —0.60 ± 0.13 ± 0.14 and 
—0.44 ± 0.51 ± 0.31 for J /ijj —> pp and J/ip — > AA, respectively. This gives an a value for 
J/i/j — > pp of 0.43 ± 0.13 ± 0.14, which is smaller than, but still consistent with, the current 
world average 0.66 ± 0.05. We also report the first observation of rj c — > AA decays with 
B(rj c -> AA) = (0.8718;^;™ ± 0.27) x 10~ 3 . The observed ratio B(r] c -> AA)/B(r) c -> pp) 
is 0.67lo;}g ± 0.12, which is consistent with theoretical expectation [lril ]. 
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